Spin Current in Spin-Orbit Coupling Systems 
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We present a simple and pedagogical derivation of the spin current as the linear response to an 
external electric field for both Rashba and Luttinger spin-orbital coupling Hamiltonians. Except for 
the adiabatic approximation, our derivation is exact to the linear order of the electric field for both 
models. The spin current is a direct result of the difference in occupation levels between different 
bands. Moreover, we show a general topological spin current can be defined for a broad class of 
spin-orbit coupling systems. 
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Spintronics aims to manipulate spins of particles. As 
such, an essential step in the field is the generation of 
a reliable spin current. The injection of spin polarized 
electron current from a ferromangetic metal is not favor- 
able because polarization is lost at the interface due to 
the conductance mismatch 0,0- Injections from ferro- 
magnetic semiconductors into nonmagnetic semiconduc- 
tors have been successfully developed in recent several 
years 0, 0, • The theory of spin transport in the fore- 
mentioned cases depends on the detailed mechanism of 
spin relaxation, where spin transport generally is a dissi- 
pative process. 

Recently, Murakami, Nagaosa and ZhangQ have dis- 
covered a dissipationless and topological spin Hall current 
in the hole doped semiconductors with strong spin orbit 
coupling. These authors studied the effective Luttinger 
Hamiltonian 0] 
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which describes conventional semiconductors such as Si, 
Ge, GaAs and InSb. The geometrical structure of the 
effect is such that for an electric field applied on the z 
direction, a y-polarized spin current will flow in the x 
direction. The electric field induced spin current can be 
summarized by the following formula 



(2) 



where e^fc is the antisymmetric tensor. Unlike the ordi- 
nary Ohm's law, this equation has the remarkable prop- 
erty of time reversal symmetry 0, 0- This effect also has 
a deep topological origin[(|0, related to the topolo gica l 
structure of the four dimensional quantum Hall effect [13|. 

In the other independent work by Sinova et alQ, the 
dissipationless, or the intrinsic spin current is also dis- 
covered to exist in the two dimensional Rashba system, 
described by the Hamiltonian^^, 0] 
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In this case, the spin current is polarized in the direction 
perpendicular to the two dimensional plane and flowing 
in a planar direction perpendicular to the direction of 
the charge current. Surprisingly, the spin conductance 
in these systems turns out to be independent of the spin 
orbit coupling and given by: 



e 

8tt' 



(4) 



' ^{.PxSy PyS x ) 



(3) 



In both of the above models, the intrinsic spin cur- 
rent is induced by an external electric field. The authors 
°f 0- @ have argued that the spin current is dissipa- 
tionless because the spin conductance is invariant under 
time-reversal operations. However, the approaches taken 
by the authors of || and |(| are markedly different. In 
the first two-dimensional model a semi-classical ap- 
proach is used to derive the spin current, while in the 
second three-dimensional model the spin current is 
derived as a topological effect in momentum space. It 
therefore appears that the physics of the spin current in 
these two models is different. In fact, in the first model, 
there is no known topological structure. However, from a 
purely theoretical point of view, we should be able to de- 
rive the results and understand the physics in one unified 
approach. 

The aim of this article is to present a simple, ped- 
agogical and unified derivation of the spin current for 
both models. The derivation of Murakami et al. 0, @ 
emphasizes on the momentum space topology which re- 
quires some advanced mathematical knowledge. The 
semi-classical derivations presented in 0, 0, 0] may not 
be familiar to general readers. In view of the impor- 
tance of the effect, we feel that the general reader could 
benefit from a simple derivation based on the standard 
time dependent perturbation theory within single parti- 
cle quantum mechanics. As in |6|, |8| , we work within the 
adiabatic approximation. Since there is no interaction 
between particles, the result derived here is exact as far 
as only the linear response to the external electric field in 
the adiabatic limit is concerned. Moreover, general con- 
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elusions about the spin current can be manifestly drawn 
in this approach. Following the common definition of the 
spin current, we show that the spin current is always a di- 
rect result of the difference in occupation levels between 
different bands in the models. Part of the spin current 
can be interpreted as topological spin current. In fact, we 
present an analysis and derivation of the topological spin 
current for a broad class of spin-orbital coupling models. 

In the following, we first calculate the spin current in 
the Rashba and Luttinger models and then discuss the 
topological part of the spin current. For the Luttinger 
Hamiltonian, we show that the expectation value for the 
spin current in the heavy hole and light hole states differs 
by exactly a minus sign. This result leads to the conclu- 
sion that the contributions to the spin current from the 
heavy hole and light hole bands should be exactly op- 
posite and differ by the Fermi velocities of two bands. 
The total spin current gives a quantum correction to the 
semiclassical result of Ref. Q. The nature of the quan- 
tum correction can be manifestly understood in our cal- 
culation, and has also been discussed in Ref. Q . For the 
Rashba Hamiltonian our approach gives the same spin 
current as the Ref. ,8j if the same definition of the spin 
current is taken. 

In our calculation, the spin current is defined by the 
velocity times the spin, which is a rank two tensor. How- 



ever, the velocity operator in general does not commute 
with the spin operator in a model with spin-orbit cou- 
pling. In order to define the spin current tensor as a 
Hermitian operator, we have to symmetrize it: 
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Let us consider a general spin-orbit coupling model de- 
scribed by the many-body Hamiltonian H(P,S). In the 
presence of a constant external electric field, we choose 
the vector potential A — —Et. The total Hamiltonian 
becomes time dependent, H(t) — H(P — eEt,S). Let 
\G,t > be an instantaneous ground state eigenstate of 
the time-dependent Hamiltonian, 



H(t)\G,t >=E G (t)\G,t > 



(6) 



The many-body ground state wavefunction ^c(t) of the 
Hamiltonian satisfies the Schrodinger equation, 



i^a(t) = H(t)V G (t). 



(7) 



By first-order time-dependent perturbation theory, we 
have 



|* g (t)) = exp{-* jf dt'E G (t')}{\G, t) + < £ j (1 - e W«)-*>ffl)«)}, 
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where \n, i) are the instantaneous excited eigenstates. 

The second term in the RHS of the expression above contains a fast oscillation term which averages to zero, and 
which we neglect below. For the non-interacting Fermi system, the above expression can be simplified into summation 
over all instantaneous single particle eigenstates. Then for an arbitrary operator O, the difference of expectation 
values between the perturbed and unperturbed states is given in [l^j 

/lT , „a UT( A . sr {x,P(t)\d\\',p(t))(\',p(t)\§- t \\,p(t)) + (\,p(t)\jL \x',p(t))(x',p(t)\d\x,p(t)) 

« sws, e \' ,P{t) - e \,P{t) 



where |A, Pit)) is the instantaneous eigenstate with po- 
larization A of the single particle Hamiltonian. 



The entire calculation of the spin current that follows 
relies on the above Kubo formula. However, when the 
Hamiltonian has degenerate states, we use the following 
convention: if a set of states, {|A, P(t) >}, are degenerate 
in energy, we can always choose a complete orthogonal 
basis of states in the set, {\a,P(t) >}, such that, for 
any two new different orthogonal states \ctx,P(t) > and 



\ot2,P(t) >, we have, 

( ai ,P(t)\^\a 2 ,P(t)) = 0. (10) 

In this case, the summation index in the formula does not 
include the degenerate states and therefore, the formula 
is well defined. 

Let us now consider the particular case of the Rashba 
Hamiltonian in an external electric field. The time de- 
pendent Rashba Hamiltonian is given by 

P(t) 2 

H R {t) = + l{P x {t)S y - P y (t)S x ). (11) 
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For a given P(t), the instantaneous eigenstates are given 
by 



\X,P(t) >=U R \X>, 



(12) 



with e(P(t)) = ^ + lX\P(t)\, U R = e 



,0W 



tan 



P«(i) 



the azimuthal angle, and |A) is the eigenstate 



of S y with S y \X) = A|A). Therefore, 



||A,P W > 



-f^5,|A,P(t)>, (13) 



with = ee ijEi p 3 -i(t) ) where is rank- 2 antisym- 

metric tensor. By applying the Kubo formula, for any 
Hcrmitian operator O, we obtain 

p. (7 s ) 

(A,P(f)|0|A,P(i)>=2ee^-^ 
x E (^P W |^|A,P(t))MM^P^. (14) 
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For the spin one half particles, the above formula is sim- 
plified to 
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±-,P(t)\0\ ± -,P{t)) = ±C a e t3 E, 



p{tf 



(15) 



where C = e/^Re(^\0(t)\ — |). The spin current op- 
erator when the spin is polarized in the perpendicular 
direction to the xy plane is given by 



Pi 7 
Ji = — S z + —eij{S z Sj + SjS z ). 
m 2 



(16) 



It is easy to show that Pe(±|Jj(i)| - 1) = ^(H^l-f) = 
Considering the whole fermi surface, we can easily 
calculate the total spin current. Let's take the electric 
field is in x direction and consider the current J y . 
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where the APf is the difference of the fermi velocity for 
the two bands. In this model, APf = mX. This yields 
the result of the EqQ] 

We now turn to the discussion of the effective Luttinger 
Hamiltonian. In the presence of the external electric field, 
the time dependent effective Luttinger Hamiltonian is 

H L (t) = i(7i + \l2){P{t) f + 2 72 (P(t) • S) 2 }. (18) 
2m 2 

For a given P(t), the Hamiltonian has four instantaneous 
eigenstates, 



where e x (P(t)) = ^f*( 7l + (f - 2A 2 ) 72 ). For A = ±| 
and A = ±i, they are referred to as the heavy hole band 
and light hole band respectively. The eigenstates can be 
explicitly written as 

|P(t), A >= U L \X >, U L = e -^(*)^ e - ie (*) s «|A >, (20) 

where tan(0(t)) = P y {t)/P x {t), cos(<9(*)) = P z (t)/\P(t)\ 
and S Z \X >= A|A >. 

Since the eigenstates are degenerate, we have to choose 
an orthogonal basis and satisfy the Eq.^|in order to use 
the Kubo formula. Without loss of generality, we choose 
the electric field in the z direction. In this case, <fi is 
time-independent. Therefore, 

(P(t),X'\^\P(t),X)-i^.(X'\S v \X) (21) 

From this equation, we obtain that for the states 
with the hclicity equal to ±|, the matrix element, 
(P(t) 7 — ||9t|P(t), |) vanishes. The only states for which 
we have to find an orthogonal base are in the helicity ± i . 
Let us define 

\P(t),+) = ^=\P(t),^+i\P(t),~), 

\P(t),-) = ^={\P{t),\)-i\P{t)-\), (22) 

which satisfy (P(t), +\d t \P(t),~) = 0. 

For an arbitrary operator O, let 0(t) = U^OUl- By 
applying the Kubo formula, we obtain the expectation 
value for an arbitrary Hermitian operator, 

(*,(*)|d|*|(t)) = ^^Im«i|d(t)|§» (23) 



and 



(* + (t)|0|*+(t)) = 



V3m d6(t) 



W « {M<ll6(*)|i> + (-i|6(«)|-|) 

+Re({-I|6(t)|) + (i|d(t)|-|))} (24 



The spin current operator where the spin is polarized in 
y direction and flows in x direction is given by 



7i + §72 



PxSy - ¥-[S y ((P ■ S)S X + S X (P ■ S)) + h.cp5) 



Jy 2m ^ 

The matrix element Im(|| J x {t)\\) is calculated to be 

Im<|l^(*)l^) - ^ sin0( 7l cos 2 + 2 72 sin 2 0). (26) 
We thus obtain, 



P(t)|P(t),A>=e A (P(t))|P(*),A>, 
P{t)-S, 



\P(t)\ 



-\P(t),X>=X\P(t),X). 



(19) 



< *»(*)|J*|*3(*) >= - < >= 
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(71 Pi + 2 72 P V 2 )P 



(27) 
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where we used d0(t)/dt = eEsm(0(t))/\P(t). 

For the states \P(t), — |) and \P(t),—), the expecta- 
tion values for spin current are the same as \P(t), |) and 
\P(t),+) respectively. Therefore the total spin current 
by including all the particles in the two bands is given by 



n = eig( ? 2 + 272) AP/, 

4tt 2 72 1 



(28) 



where APf is the fermi momentum difference between 
the heavy and light hole bands. Once again, we show 
that the spin current comes from the occupation differ- 
ence between two bands. The above result has been in- 
dependently obtained by Murakami, Nagaosa and Zhang- 
through a slightly different derivation based on Kubo for- 
mula too0- 

This full quantum mechanical calculation gives a quan- 
tum correction to the original semiclassical result 6] . In 
fact, the difference comes from the definition of spin cur- 
rent operator. In |6(], an effective Hamiltonian was de- 
rived by introducing a monopole in momentum space. 
The spin current is thought of as a topological effect of 
the monopole. In the heavy hole states, the gauge po- 
tential is abelian while in the light hole states, the gauge 
potential is nonabelian. However, the field strength in 
both bands is abelian. For each helicity states, the field 
strength is given by 



[A, Dj] = A(A 2 
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(29) 



where the operator Dj = P h ^(U^d Pj UL) (A(2) is the 
projection onto heavy (light) hole bands). This gauge 
field modifies the semiclassical equation of motions as 



Pi 

Vi,x = Xt \ = h F ik E k . 



(30) 



In the spin current is derived by replacing the spin 
operator by its expectation value in the helicity states 
compared to the definition in the Eq[S] If we use the 
same replacement, it is straightforward to show that the 
spin current from the perturbation theory is the same as 
from the 0- Namely, for a given helicity state |A, P(t) >, 
the expectation value of < Uj >\ in the adiabatic ap- 
proximation from the perturbation theory is given by 
^<A|^|A>. 

However, the above calculation does not underline the 
topological nature. Several questions still remain. The 
first is how the topological spin current can be separated 
from the general spin current formula in the Eq0 The 
second is that since the calculation is performed on a 
specific model, it is not clear whether the topological ar- 
guments can be applied to more general cases such as 
realistic anisotropic, inversion-symmetry breaking semi- 
conductors. The recent work of Murakami Nagaosa and 
ZhangjTj have answered some of these questions from the 
Kubo formula for the isotropic Luttinger Hamiltonian. 



Here we give an independent argument based on our for- 
malism. 

Let us review Eq^land discuss a general case. Let's as- 
sume a general unitary transformation U which is a func- 
tion in the momentum space and diagonalizes a general 
Hamiltonian H . For any operator O, let 0(U) — U + OU . 
Therefore, Hq = H(U) is diagonal and 



J!(U) = Y [S i (U)[X j (U),H ] 



[X J (U),H }S l (U)}(31) 



Let's write S l (U) = Sf(U) + Sf(U), and X 3 (U) = 
X?(U)+X?(U) where O p (U) keeps the elements of 0{U) 
which are only between the degenerate eigenvalues of 
H(U) for a given operator O. Namely it is the projec- 
tion onto the degenerate bands. O c (U) is the leftover 
part. Now we can define the total spin current operator 
into J{{U) = T?(U)+A{(U), where the first part Tf (U) 
is defined 

niU) = ^([SfWXjW+XjSfiUlHo] 

+ [Si(U)X?(U)+X?Si(U),H ]) (32) 



and 



Ai(U) = Y [S?(U)[X°(U),H ] 



[x^u),H }smm 



where the relations [Sf(U), H ] = and [X?(U),H ] = 
have been used in the above equations. It is clear that 
A\ ({/) is the band crossing contribution to the spin cur- 
rent. T-(U) can be considered as the topological part of 
the spin current. This statement is true for any mod- 
els with arbitrary number of bands and with arbitrary 
degeneracy in each bands caused by spin orbit coupling. 
The proof is straightforward from the perturbation the- 
ory. 

Without loss of generality, we assume that Hq is the 
( EJ mi ...\ 
diagonal matrix, Ho = £ , 2-^m 2 ••• j where 

mi, m2, ■■■ are the number of degeneracies of each of the 
bands. By a direct calculation from the Eq[5] , the spin 
current contribution from T? (U) is given by 



< T( >= zTr{Sf(U)[XnU),X p k (U)}} 



dP k (t) 
dt ' 



(34) 



The above equation is independent of U for all of unitary 
matrixes which U + HU = Hq. Therefore, the symmetry 
group for U is SU(mi) SU{m 2 ) .... The above for- 
mula is manifestly gauge invariant if we view the sym- 
metry group as a gauge group in momentum space as 
described in 0. 

From the above analysis, we see that the topological 
spin current exists in much broader spin-orbit coupling 
systems. However, it requires the degeneracy of the 
bands, namely the non-abelian gauge structure in mo- 
mentum space. The direct consequences from this result 



5 



is that for a realistic anisotropic Luttinger Hamiltonian, 
the topological part of spin current will still exist, and 
that for the Rashba Hamiltonian there is no topological 
part of the spin current since there is only a U(l) gauge. 
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